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ABSTRACT 

We review the theory of electron-conduction opacity, a fundamental ingredient in the computation 
of low-mass stellar models; shortcomings and limitations of the existing calculations used in stellar 
evolution are discussed. We then present new determinations of the electron-conduction opacity in 
stellar conditions for an arbitrary chemical composition, that improve over previous works and, most 
importantly, cover the whole parameter space relevant to stellar evolution models (i.e., both the regime 
of partial and high electron degeneracy). A detailed comparison with the currently used tabulations 
is also performed. The impact of our new opacities on the evolution of low-mass stars is assessed 
by computing stellar models along both the H- and He-burning evolutionary phases, as well as Main 
Sequence models of very low-mass stars and white dwarf cooling tracks. 

Subject headings: conduction - stars: horizontal-branch - white dwarfs - stars: interiors - stars: 
evolution 
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1. INTRODUCTION 

One of the main physical inputs required to solve the 
equations of stellar structure is the electron thermal 
conductivity. When the degree of electron degeneracy 
is significant, electron conduction is the dominant en- 
ergy transport mechanism, and the value of the electron- 
conduction opacity - proportional to the inverse of the 
electron thermal conductivity - enters the equation of the 
temperature gradient. This physical condition is verified 
in the interiors of brown dwarfs, very low-mass stars with 
mass Mtot < 0.1 5M© (the exact value slightly affected by 
the metal content; see Chabrier & Baraffe 2000 for a re- 
view) , in the He-core of low- mass stars during their Red 
Giant Branch (RGB) evolution (see, e.g., Salaris, Cas- 
sisi, & Weiss 2002), in the CO core of Asymptotic Giant 
Branch stars, as well as in white dwarfs (inside the core 
and in a portion of the envelope - see, e.g., Prada Moroni 
& Straniero 2002) and envelopes of neutron stars. 

Until the mid-1990s, the main sources of electron- 
conduction opacity for stellar model computations were 
those provided by Hubbard & Lampe (1969), by Itoh and 
coworkers (Flowers & Itoh 1976, 1979, 1981; Itoh et al. 
1983, 1984; Mitake, Ichimaru, & Itoh 1984; Itoh & Ko- 
hyama 1993; Itoh, Hayashi, & Kohyama 1993), and by 
Yakovlev and coworkers (Yakovlev & Urpin 1980; Urpin 
& Yakovlev 1980; Raikh & Yakovlev 1982; Yakovlev 1987; 
Baiko & Yakovlev 1995). Each of these sources of opac- 
ities has its own shortcomings. The Hubbard & Lampe 
tabulations cover a very limited set of element mixtures, 
do not take into account relativistic effects, and contain 
significant gaps in the temperature-density diagram (see 
§2.3 below). The formulae by Itoh et al. (1983, 1984, 
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1993) and Itoh & Kohyama, as well as Yakovlev & Urpin, 
Raikh & Yakovlev, and Baiko & Yakovlev, are succes- 
sive improvements - in their domains of validity - over 
the Hubbard & Lampe results for the opacity contribu- 
tion due to electron-ion (ei) scattering: they take into 
account relativistic effects, more accurate structure fac- 
tors for an ion liquid, phonon-electron interactions for an 
ion solid, and can be employed to compute conductive 
opacities for arbitrary chemical mixtures. Simple formu- 
lae for the electron-electron (ee) scattering contribution 
to the conductive opacities have been derived by Urpin 
& Yakovlev (1980) and Potekhin, Chabrier, & Yakovlev 
(1997). However, the above-cited results by Itoh's and 
Yakovlev's groups are applicable only when electrons are 
strongly degenerate (T <C Tp, where Tp is the Fermi 
temperature) a condition which is not strictly satisfied 
in RGB stars (see, e.g., Fig. 18 in Catelan 2005). In 
addition, they do not take into account important col- 
lective plasma effects near the solid-liquid phase transi- 
tion (Baiko et al. 1998). These effects have been taken 
into account by Potekhin et al. (1999), whose calcula- 
tions can be also used to compute opacities for arbitrary 
astrophysical mixtures at T < Tp. For the ei scattering 
contribution, the latter calculations have been extended 
to partially degenerate plasmas by Ventura & Potekhin 
(2001), using a thermal averaging procedure presented 
by Potekhin (1999) (see § 2.3 below). 

The value of the electron-conduction opacity is crucial 
in particular for the evolution of low-mass RGB stars and 
the fate of their progeny (see, e.g., the extensive discus- 
sions in Catelan 2005 and Catelan, de Freitas Pacheco, 
& Horvath 1996). The whole thermal stratification in- 
side the He-core and the temperature at the base of the 
surrounding H-burning shell depend strongly on the ef- 
ficiency of electron conduction (for a detailed discussion 
on this issue see Catelan et al. 1996; Salaris et al. 2002, 
and references therein). As a consequence, the mass size 
of the He core (M c h c ) at the He-burning ignition (the 
He-flash at the RGB tip) depends on the value of the 
electron conductivity. An accurate determination of the 
M C He value is extremely important for several astrophysi- 
cal problems. First of all, the value of M c r- c regulates the 
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stellar luminosity at the RGB tip and, in turn, the evo- 
lutionary lifetime along the RGB: any increase of M c h c 
causes an increase in both the surface luminosity at the 
RGB tip and the RGB evolutionary lifetime. Given that 
the mass loss efficiency along the RGB increases with in- 
creasing surface luminosity and decreasing effective tem- 
perature (see, e.g., Reimers 1975) the amount of mass 
lost during the RGB evolution is strongly influenced by 
any variation in the mass of the He core at the He-flash. 
This has a sizeable effect on the morphology of the Hor- 
izontal Branch (HB) (see the review by Renzini & Fusi 
Pecci 1988 for a detailed discussion on this issue). In 
addition, one has to notice that the brightness of the 
RGB tip is a commonly used standard candle for old, 
metal-poor stellar systems (Lee, Freedman, & Madore 
1993). Changes in the theoretical predictions of M cHe 
have a direct influence on the theoretical calibration of 
this distance indicator and, in turn, on the derived dis- 
tances (Salaris & Cassisi 1998). Also, the luminosity of 
the core He-burning phase in low-mass stars is mainly 
controlled by the value of M c h c - Given that the HB lu- 
minosity and, in particular, the mean brightness of RR 
Lyrae stars is one of the most important distance indi- 
cators for old stellar systems, the theoretical calibration 
of this standard candle depends on the accuracy of the 
electron-conduction opacities. 

The value of M c h c also affects the duration of the core 
He-burning phase, playing therefore a role in the theo- 
retical calibration of the so-called i?-parameter (number 
ratio of HB stars to RGB stars brighter than the HB), 
which is an indicator of the initial He content of old stel- 
lar populations (Iben 1968). When using the Itoh et al. 
conductive opacities, Castellani & Degl'Innocenti (1999) 
found an increase by ~ 0.005 M Q of M c n e at the He flash 
for a 0.8 M Q star with metallicity Z = 2 x 10 -4 compared 
to models based on the Hubbard & Lampc (1969) opac- 
ities. The same numerical experiment performed on a 
1.5 Mq star with initial solar composition showed an in- 
crease of ~ 0.008 M©. A similar comparison performed 
with the Potekhin et al. (1999) opacities provides He- 
core masses at the He-ignition between those obtained 
with the Itoh et al. and the Hubbard & Lampe ones, but 
closer to the Itoh et al. results. 

The above-cited papers by Itoh et al., Yakovlev et al., 
and Potekhin et al. were devoted mostly to plasmas at 
the conditions typical for neutron star envelopes, where 
ee scattering is relatively unimportant. Therefore the ee 
scattering mechanism was paid little attention in those 
papers. For instance, Potekhin et al. (1997, 1999) used 
a fitting formula for the ee contribution, valid only for 
strongly degenerate electrons. However, as already noted 
by Hubbard & Lampe (1969) and recently stressed by 
Catclan (2005), the ee scattering can give a consider- 
able contribution to the opacity in partially degenerate 
regions within the He-core of RGB stars. In this paper 
we improve the treatment by Potekhin et al. (1999) and 
Potekhin (1999) for the case of a non-magnetized plasma 
by including the ee scattering in partially degenerate and 
nondegenerate matter. In addition, we take into account 
a recent improvement of the treatment of the ee scatter- 
ing in degenerate matter by Shternin & Yakovlev (2006). 
To verify the impact of these new conductive opacities, 
we compute models of low-mass stars and perform a de- 
tailed analysis of the effects of the new conductive opac- 



ity evaluations on both the RGB and HB evolutionary 
phases. We also assess their impact on models of very 
low-mass stars and white dwarfs. 

The paper is organized as follows: in § 2 we briefly 
discuss the physical mechanisms governing the electron 
conduction in the different physical regimes and empha- 
size the assumptions made by the various existing treat- 
ments - including our new updated results - with the 
aim of highlighting the limitations of the available tabu- 
lations of conductive opacities. An analysis of the impact 
of our new conductive opacities on stellar models is pre- 
sented in § 3. A summary and final remarks follow in 
§ 4. 

2. STELLAR CONDUCTIVE OPACITIES: THE 
STATE-OF-THE-ART 

Let us recall that the total opacity of stellar matter 
(e.g., Marshak 1940; Carson 1976) can be written as 

fttotai = V( K rad + K c 1 )^ where K rac j and k c are the radia- 
tive and conductive opacities, respectively. The latter is 
related to the thermal conductivity k by the equation 
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where a is the Stcfan-Boltzmann constant, T is the tem- 
perature, and p is the density. 

The kinetic method (Ziman 1960) is most practical for 
calculating >c. Using the elementary theory in which the 
effective electron scattering rate v does not depend on 
the electron velocity, one can write (Ziman 1960) 



n e k 2 T 



x = a- 



(2) 



where n e is the electron number density, m e is the elec- 
tron mass, k is the Boltzmann constant, a = 3/2 for a 
nondegenerate electron gas (T ^> Tp), and a = 7r 2 /3 for 
strongly degenerate electrons (T <C Tp). If the electrons 
are degenerate and relativistic, one should replace m e in 
Eq. (2) by m* = m e 7 r , where 7 r = \J\ + x 2 , 



= p F /m e c = 0M009{pZi/A) 



1/3 



(3) 



is the relativistic density parameter, pp = S(37r 2 n e ) 1 / 3 is 
the Fermi momentum, Zi and A are the nuclear charge 
and mass numbers, and the mass density p is expressed 
in g cm" 3 . Note that the electron Fermi temperature 
that accounts for relativistic effects can be expressed as 

Tp = (to* - m e )c 2 /k = 5.93 x 10 9 K ( 7r - 1) , (4) 

where k is the Boltzmann constant. 

In a fully ionized plasma, v is determined by electron- 
ion and electron-electron Coulomb collisions. Let us as- 
sume that the effective frequencies of different kinds of 
collisions simply add up, i.e., v = v e i + v ee . This so-called 
Matthiessen rule is strictly valid only for extremely de- 
generate electrons (see Hubbard & Lampe 1969). How- 
ever, in practice it gives a good estimate of the conduc- 
tivity: one can show that v e i + v ee < v < v e i + v ee + 8u, 
where 5v <C rnin^ej, z/ ee ) (sec Ziman 1960). 

2.1. Nondegenerate Electron Gas 

Let us consider first the case where the electrons are 
non-degenerate (T ^> Tp) and non- relativistic (x r <C 1). 
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This case has been studied, e.g., by Spitzcr & Harm 
(1953), Braginskii (1958), and Spitzer (1962). The ef- 
fective energy-averaged ei collision frequency is 

where rii is the ion number density, mi is the ion mass, 
and A e j is the so-called Coulomb logarithm. In the con- 
sidered case A e j is a slowly varying function of density 
and temperature. In general, its value depends on the ap- 
proximations used to solve the Boltzmann equation, but 
in any case its order-of-magnitude estimate is given by 
the elementary theory, where the Coulomb collision in- 
tegral is truncated at small and large impact parameters 
of the electrons. Then one obtains A ei ~ ln(r max /r m ; n ), 
where r max and r m j n are the maximum and minimum 
electron impact parameters. The parameter r max can be 
set equal to the Debye screening length, r~ ax = 47r(n e + 
Zfm)e 2 /kT. The second parameter can be estimated 
as r m ; n = max(AT, r c \), where At = jmJtT is 

the thermal de Broglie wavelength, which limits r m i n in 
the high-temperature regime (where the Born approxi- 
mation holds), and r c \ = Zie 2 /kT is the classical closest- 
approach distance of a thermal electron, which limits 
r m ; n in the low-temperature, quasiclassical regime. 
A similar effective frequency 

characterizes the efficiency of ee collisions. If A ee ~ A e j, 
then v e i/v ee ~ Zj, therefore for large Zi the ei collisions 
are much more efficient than the ee collisions. 

2.2. Strongly Degenerate Electrons 

The thermal conductivity of strongly degenerate elec- 
trons in a fully ionized plasma is given by Eq. (2) with 
a = 7r 2 /3. In order to determine the effective collision 
frequency v that enters this equation, we shall use the 
Matthiessen rule, mentioned above. This will allow us to 
consider the different types of electron scatterings sepa- 
rately in the following sections, and calculate v as their 
sum. 

2.2.1. Electron-Ion Scattering 

Thermal transport coefficients of degenerate electrons 
were studied in a number of papers. In the 1990's, the 
formulae for conductive opacities mostly used in astro- 
physics were those derived by Itoh and coworkers and 
by Yakovlev and coworkers (see references in § 1). These 
formulae are valid for degenerate electrons in strongly 
coupled ion liquids or solids - i.e., at T -C Tp and T > 1, 
where 

(Z,e) 2 _ 0.2275 2 / P y/3 

is the Coulomb coupling parameter 5 (here a, = 
[ATmjS]" 1 ^ is the ion sphere radius, T 6 = T/10 6 K, 
and p is expressed in g cm -3 ). 

5 This parameter regulates the state of the Coulomb plasma: 
the ions constitute a gas at T < 1, a liquid at 1 < T < r m , and a 
crystal at V > r m , where r m fs 175 (see, e.g., Potekhin & Chabricr 
2000). 



It is convenient to present the effective frequency of ei 
collisions in the form 

Vei = 5 UiAei ( 8 ) 

Pp1>F 

(where vf — pp/ml = cx r /~f r is the electron Fermi ve- 
locity), and derive a fitting formula for the Coulomb log- 
arithm A e j function of n e , T, and Z{. Yakovlev 
& Urpin (1980) used Eq. (8) in the ion liquid regime, 
where A ei is a slow function of order unity. They em- 
ployed another form of v ei for an ion crystal, where col- 
lective effects are important - that is, where electrons are 
scattered not by individual ions, but by their collective 
excitations, the phonons. 

However, Baiko et al. (1998) showed that in fully ion- 
ized dense plasmas there is no appreciable discontinuity 
between the conductivities in the solid and liquid phases. 
This is due to the incipient long-range ion structures in 
the ion liquid and multiphonon scattering in the solid 
near the melting point Y ~ T m . Thus, the conductivity 
is smooth across the melting. More recent and more ac- 
curate calculations (Potekhin et al. 1999) resulted in new 
fitting formulae for the ei conductive opacities in degen- 
erate matter. These calculations included the accurate 
numerical structure factor for a Coulomb liquid calcu- 
lated by F. Rogers & H. E. DeWitt (unpublished) and 
fitted by Young, Corey, & DeWitt (1991), for a Coulomb 
coupling parameter T > 1. Potekhin et al. (1999) have 
taken account of the modifications of the effective (for 
electron scattering) ion structure factor around the melt- 
ing point, discussed by Baiko et al., as well as the quan- 
tization of ion motion at T below the ion plasma temper- 
ature T Pj ; on and "freezing-out" of the umklapp processes 
at T < T pAon (Raikh & Yakovlev 1982). A correct lim- 
iting form of the ion structure factor at T <C 1 was also 
ensured. Hence the results of Potekhin et al. (1999) are 
probably valid at any T. On the other hand, it should 
be noted that rigorous calculations for T < 1 have never 
been carried out, and would certainly prove of interest in 
placing this "smooth interpolation argument" on a more 
solid physical basis. 

The absence of a big discontinuity at T — r m allowed 
Potekhin et al. (1999) to construct a single fitting formula 
for A ei using Eq. (8) at all temperatures. We shall follow 
this approach hereafter. It should be noted, however, 
that A e i in Eq. (8) is a slowly varying function only at 
T < r m , but it rapidly decreases in the solid regime at 
T > T m (see Potekhin et al. 1999 for details). 

2.2.2. Electron- Electron Scattering 

The expression of u ee for the relativistic degenerate 
electrons at T <C T p was obtained by Flowers & Itoh 
(1976). Here T p — huj p /k is the electron plasma tem- 
perature, determined by the electron plasma frequency 
uj p = y/<iire 2 n e /m*. Urpin & Yakovlev (1980) extended 
the results of Flowers & Itoh (1976) to higher temper- 
atures, where T C Tf, but not necessarily T <C T p . 
Potekhin et al. (1997, 1999) calculated ee conductive 
opacities according to the theory of Urpin & Yakovlev 
(1980) and presented a fitting formula for their results. 

Recently, Shternin & Yakovlev (2006) have reconsid- 
ered the problem taking into account the Landau damp- 
ing of transverse plasmons. This effect is due to the 
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difference of the components of the polarizability tensor, 
responsible for screening the charge-charge and current- 
current interactions: the transverse current-current in- 
teractions undergo "dynamical screening." It was ne- 
glected by Urpin & Yakovlev (1980) but later studied by 
Heiselberg & Pethick (1993) in the context of the trans- 
port coefficients of the quark-gluon plasma. Shternin & 
Yakovlev showed that the Landau damping of transverse 
plasmons strongly increases v K 
and T <C T p . Their result can be written as 



in the domain of x r > 1 



3/2 



m e c 6a f 
Vee = — — ^jj- x r y VA- I(/3 T ,y), 

where at is the fine-structure constant, /3 T ■ 
x r / 7r , y = V3T p /T = (571.6/T 6 ) v / /5;x r , and 

8/315 \ 
+ 0.0435y J 



(9) 



vf/c 



7 ^) = Ki 



+0 



+ 



x In | 1 + 
, 2.404 



128.56 



B 



37.1y + 10.83y 2 + y 3 
C- 2.404/5^ 



x In 



l + D 



1 + 



C 



1 + 0.1/32/ 
B 



x In 



1 



A(3y + {Pyf 
18.52/? 2 L> 
B 

B 



Ay + 10.^{(iyY + {f3yf/\ ' 

12.2 + 25.2 /? 3 , B = Acxp{(0.123636 + 
C = 8/105 + 0.05714 f3 l , and D = 



where A = 
0.016234 /3 2 )/C], 
0.1558 y^-^P. 

These equations are valid only for T < Tp; their ex- 
tension to the case T > Tf will be suggested in § 2.3. 

It is easy to check that Eqs. (9) and (10) reproduce 
the low-density (x r c 1) asymptote of the previous fit to 
v ee presented by Potekhin et al. (1997, 1999). However, 
the deviation from the previous results is already large 
at x T <~ 1 and can reach many orders of magnitude at 
x r > 1. 

2.2.3. Impurities and Mixtures 

Real stellar plasmas are often mixtures of different 
chemical elements. In this case Eq. (3) for the relativity 
parameter x r has to be modified by replacing Zi/A by 
the unweighted mean value of the charge-to-mass ratio, 
(Zi/A), averaged over all species. The effective collision 
frequency v e i should also be modified. The required mod- 
ification can be different, depending on the state of the 
plasma and on the amount of impurities. For example, 
Flowers & Itoh (1976), Yakovlev & Urpin (1980), and 
Itoh & Kohyama (1993) considered electron scattering by 
charged impurities in a Coulomb crystal. If the fraction 
of impurities is small and they are randomly distributed, 
electron-impurity scattering can be treated as scattering 
by charge fluctuations, controlled by the impurity pa- 
rameter Q = J^j Yj (Zj ~ {^)) 2 i where Yj = nj/ J^j n j 1S 
the number fraction of impurities of the jth kind, Zj is 
their charge number, and (Z) is the mean charge number, 
which in the considered case is close to the charge num- 
ber of the main ion species that forms the crystal lattice. 



Then, using the Matthicsscn rule, one can obtain i/ e , as 
a sum of the terms corresponding to the ei (electron- 
phonon) scattering in a homogeneous lattice and to the 
electron scattering by charge fluctuations, which is given 
by Eq. (8) with Zi replaced by Z- lmp = y/Q. This ap- 
proach has been adopted by Yakovlev & Urpin (1980), 
Potekhin & Yakovlev (1996), and Gncdin, Yakovlev, & 
Potekhin (2001), and implemented in the online database 
referenced below. 

An alternative approach is relevant when there is no 
dominant ion species which forms a crystal (e.g., in a 
liquid, a gas, or a glassy alloy). In this case, one can use 
Eq. (8) with ZfriiA ei replaced by J2j ^j n j^iiy where 
summation is over all ion species j, and the Coulomb 
logarithm A J ei depends generally on j. An approximation 
to A J ei based on the plasma "additivity rule" has been 
suggested by Potekhin et al. (1999). A much simpler yet 
reasonable approximation to the conductive opacity can 
be obtained by using an effective charge number equal to 

\[W) = (Ej^S) 1/2 - 

A brief discussion and comparison of the two ap- 
proaches has been given by Brown, Bildsten, & Chang 
(2002). 

In the present paper, we consider pure He composition 
for the core of RGB stellar models. This is justified by 
low mctallicities in most of these models, which make 
the exact conductivity of a mixture close to the pure He 
conductivity. For instance, an admixture of 0.1% (by 
number) of N or C to the He liquid raises the conductive 
opacity by w 1%, which can be safely neglected. 

2.3. Partially Degenerate Electrons 

The thermal conductivity of partially degenerate elec- 
trons is more difficult to calculate. For the fully ionized 
Coulomb plasmas this problem was studied by Lampe 

(1968) , who calculated the transport coefficients of a 
weakly coupled plasma of nondegenerate ions and par- 
tially degenerate nonrelativistic electrons by means of 
a Chapman-Enskog solution of the quantum Lenard- 
Balescu kinetic equation, including both ee and ei colli- 
sions. The employed method is more fundamental and 
accurate than the simple estimates presented in Sect. 2.1. 
However, the ion correlations have not been taken into 
account, i.e., T <C 1 was assumed. Hubbard & Lampe 

(1969) combined these calculations with earlier results 
of Hubbard (1966), who considered the transport coeffi- 
cients determined by the ei collisions in a nonrelativistic 
degenerate electron gas, taking into account ion-ion cor- 
relations. Hubbard & Lampe provided conductive opac- 
ities in tabular form for various chemical compositions. 

In the domains of degenerate electrons and strongly 
correlated ions, more accurate results have been obtained 
after the work of Hubbard & Lampe (1969), as discussed 
in Sect. 2.2. Although the Hubbard & Lampe results 
are sufficiently accurate at T » Tp, their tabular form 
of presentation does not allow their extension to other 
chemical compositions or plasma parameters beyond the 
tabulated range. 

Potekhin (1999) proposed a procedure of averaging the 
expressions for ei transport coefficients, obtained in the 
limit T/T-p — > 0, over the thermal distribution of the 
electron energies. Although this procedure was originally 
devised for taking account of thermal broadening of os- 
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Fig. 1. — Conductive opacity (k = k c = K e ; + K ee ) of helium 
as a function of the degeneracy parameter T/Tp: comparison of 
Hubbard & Lampe (1969) data (crosses) with the present ei + ee 
results (solid lines) and with the ei opacities (dashed lines) for 
values of log T (K) marked near the lines. 



dilations of the transport coefficients in quantizing mag- 
netic fields at T -C Tp, it was shown to reproduce also 
the correct non-magnetic limit of x at T 3> Tp (Ventura 
& Potekhin 200f). According to this recipe, 



where 



x = k 2 T(<j 2 - of/do), 



X 



1) 



(11) 



(12) 



ei (e) m*(e) ( e * 

e is the electron energy, \ = ( e — 1^/kT, p is the electron 
chemical potential, A/"(e) is the number density of elec- 
tron states with energies smaller than e (i.e., the electron 
number density that would correspond to the Fermi en- 
ergy e at T = 0), and v e i{e) is the effective collision 
frequency of electrons with energy e. In a nonquan- 
tizing (or zero) magnetic field, M{e) — p 3 / (3ir 2 h 3 ) and 
m*(e) = \J m 2 e + (p/c) 2 , where p is the electron momen- 
tum. 

We adopt this approach to the ei opacity calculations 
at arbitrary T /Tp. 

There is no such procedure for the ee collisions. Note 
that a correct expression for v ee must reproduce the gen- 
eral dependence (6) at very high T or very low p (maybe 
up to a factor of a few, given the logarithmic accuracy 
of the estimates of A in §2.1). In particular, Eq. (6) 
implies v ee oc T~ 3 / 2 lnT at T — > oo, which agrees with 
an accurate asymptote at T Tp obtained by R. H. 
Williams & H. E. DeWitt (1968, unpublished) and re- 
produced by Eq. (16) of Lampe (1968). In contrast, the 
fit (10), being substituted into Eq. (9), produces a formal 
asymptote v ee oc T^ 1 InT at T —> oo. Clearly, Eq. (10) 
must be corrected at T > Tp. 

We construct such a correction, using the numerical 
data published by Hubbard & Lampe (1969), in such 
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Fig. 2. — Ratio of our new conductive opacities to Hubbard & 
Lampe (1969) data for helium as a function of the ion Coulomb 
coupling parameter T for 10 3 ' 8 < T(K) < 10 9 . Each line corre- 
sponds to a constant T value. Solid lines: nondegenerate matter; 
dotted lines: degenerate matter. Numbers near some of the curves 
are the corresponding values of logT (K). 



a way as to reproduce their tabulated conductivities in 
those p — T domains where their data can be trusted (see 
below). Specifically, we propose the following interpola- 
tion formula 6 for the case of arbitrary degeneracy: 



,dcg 



1 + t 2 



l+t + bt 2 v / T/Tp~ 



(13) 



where t = 25T/T F , and v^° s is given by Eqs. (9) and (10). 
The coefficient b = 135/-S/327T 7 « 0.434 corresponds to 
the ratio of the exact limiting expressions at T Tp and 
T > T F (e.g., Eqs. (16) and (17) of Lampe 1968). 

This approximation is compared to the numerical ta- 
bles of Hubbard & Lampe (1969) for helium in Figs. 1 
and 2. The agreement is good at T > Tp. In Fig. 1 
we also show the opacities calculated without taking ee 
scattering into account (the dashed lines). In this case, 
there is a noticeable difference at T > Tp. 

In Fig. 2 the solid and dotted lines correspond to the 
two parts of the Hubbard & Lampe (1969) tables: for 
weakly {T /Tp > 1.7) and strongly {T /Tp < 1.7) de- 
generate electrons, respectively. We have chosen this 
division because, for many T values, their tables con- 
sist of two disconnected parts, which merge roughly at 
T/Tp > 1.7. Figure 2 shows that there are two domains 
of strong disagreement between our opacities and Hub- 
bard & Lampe's at V > 1. These deviations take place 
near the p — T domains excluded from the Hubbard & 
Lampe tables because either the ions are strongly cou- 
pled, or their motion is quantized. In these cases the 
Hubbard & Lampe theory breaks down, therefore this 
disagreement does not invalidate Eq. (13). An addi- 
tional disagreement at T > 1 and relatively low tempera- 
tures comes from the correction by Shternin & Yakovlev 
(2006), discussed above. 

6 This represents an improvement over the formula used recently 
by Barriga-Carrasco & Potekhin (2006). 
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2.4. Summary of Recent Updates of the Conductive 
Opacities 

From the preceding discussion it is clear that the pre- 
vious conductive opacity evaluations had the following 
limitations. 

• Formulae by Itoh et al. (1983, 1984, 1993) and 
Potekhin et al. (1999) are relevant only at strong 
degeneracy, T C Tp. This condition is violated in 
the outer parts of RGB cores (Catelan 2005). 

• Itoh et al. (1983) relied on the ion structure factor 
calculations at T > 1. The condition T > 1 is also 
violated in some parts of RGB stars (Catelan et 
al. 1996; Catelan 2005). However, Potekhin et al. 
(1999) used more accurate evaluations of the ion 
structure factor at T > 1 and an interpolation at 
r < 1, which allowed them to consider conductivi- 
ties at arbitrary T in spite of the fact that detailed 
calculations for T < 1 have not yet been carried 
out. 

In addition, Potekhin et al. (1999) have taken ac- 
count of multiphonon scattering and incipient long- 
range order around the melting line (r = r m ) 
on the plasma phase diagram (Baiko et al. 1998), 
which led to a considerable correction to the opac- 
ities at r ~ 100 — 300, compared to Itoh et al. 
results. 

• Potekhin (1999) employed a specific thermal aver- 
aging for the contribution of the ei scattering to 
the conductive opacity (§2.3), implemented it in 
subsequent papers (e.g., Ventura & Potekhin 2001; 
Gnedin et al. 2001), and created an online database 
for computing the updated electrical and thermal 
plasma conductivities. 7 This allows one to com- 
pute conductivities determined by ei scattering at 
arbitrary degeneracy. 

However, the absence of a similar thermal averag- 
ing procedure for the ee scattering prevented one 
from obtaining accurate results in the case where 
the ee scattering can be important. In the con- 
ductivity database, two ways were employed previ- 
ously to extend the non-magnetic opacities beyond 
the domain of T <C Tp: (A) to ignore ee scattering 
altogether, or (B) to continue using the theory by 
Urpin & Yakovlev (1980). Model A was employed 
in the table of opacities and model B in the For- 
tran program presented in the database (the latter 
model was used in the analysis by Catelan 2005). 

• In model B the treatment of the ee scattering was 
based on the old theory of Flowers & Itoh (1976) 
and Yakovlev & Urpin (1980). As stressed in 
§ 2.2.2, the theory has been considerably improved 
recently by Shternin & Yakovlev (2006). 

In the present paper, we have done the following im- 
provements to the previous treatments of the contribu- 
tion to the conductive opacity due to the ee scattering: 

7 http: //www . ioffe.ru/astro/conduct/ 
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Fig. 3. — Ratio of the ee interaction opacity to the total (ee + ei) 
conductive opacity along the RGB core models of a O.SMq metal 
poor star discussed by Catelan (2005), as a function of the La- 
grangian mass coordinate (in solar units). For the ei contribution 
to the opacity, the fitting formulae of Potekhin et al. (1999) are 
used together with thermal averaging according to Eqs. (11) and 
(12). For the ee contribution, the solid lines correspond to Eqs. (9), 
(10), and (13), and the dashed lines to the older theory for degener- 
ate electrons (Urpin & Yakovlev 1980; Potekhin et al. 1999). Each 
curve of a series corresponds to a specific RGB stellar model, whose 
He core mass is equal to the value of the mass attained at the end 
of the curve. 



• In degenerate matter, the old theory (Urpin & 
Yakovlev 1980) and corresponding analytical for- 
mulae (Potekhin et al. 1997, 1999) are replaced by 
the new theory and analytical formulae by Shternin 
& Yakovlev (2006) (see § 2.2.2). 

• In partially degenerate and nondegenerate mat- 
ter, extrapolation of the ee opacity according to 
model B above is replaced by the interpolation for- 
mula (13), which correctly reproduces the high- 
temperature asymptote (6) and minimizes differ- 
ences with respect to numerical tables in the do- 
main of weak degeneracy (see § 2.3). 

Note that the latter interpolation does not concern the 
ei contribution to the opacity, which is calculated numer- 
ically according to Eqs. (11) and (12). 

We have updated previously developed Fortran code 
for thermal conductivities of electron-ion plasmas in ar- 
bitrary magnetic fields available online (see footnote 7) 
by including ee scattering in the particular case of zero 
magnetic field, taking into account the recent progress 
in the treatment of the ee conductive opacity, summa- 
rized above. We have also updated the table of thermal 
conductivities of fully ionized electron-ion plasmas, cal- 
culated using this code and available at the same URL. 
The table covers densities p from 10~ 6 to 10 9 g cm -3 , 
temperatures T from 10 3 to 10 9 K, and ion charge num- 
bers Zi from 1 to 60. For convenience of potential users, 
the table is supplemented by an interpolation routine in 
Fortran. 

In Fig. 3 we show the ratio of the ee interaction opac- 
ity to the total conductive opacity in the RGB core mod- 
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Fig. 4. — Ratio of the ee + ei conductive opacity k c computed 
in this paper to previous evaluations by Hubbard & Lampe («hl> 
dot-dashed lines), Itoh et al. (kt84, dashed lines), and Potekhin et 
al. («pg9, solid lines) for the same stellar models of Fig. 3. 



Fig. 5. — Differences in the He-core mass, RGB tip brightness 
and ZAHB luminosity as a function of the heavy element abun- 
dance, between stellar models computed with the 'old' conductive 
opacities by Potekhin (1999) and our new k c estimates (see text 
for details). 



els discussed by Catelan (2005). The new results (solid 
lines), which correspond to Eq. (13), are compared with 
the old model B (dashed lines, which is similar to Fig. 19 
in Catelan 2005). The differences between the new and 
old results are caused, first, by the better treatment of v ee 
in the regime of partial degeneracy according to Eq. (13), 
which decreases the opacity at T Tp but increases it 
at T ~ Tp, and secondly, by the improved treatment of 
i/fg S , reviewed in § 2.2.2, which increases the opacity at 
high densities. The opacity decrease is noticeable only 
in the outer layers of the RGB stellar models with rel- 
atively high Mtot, whereas the increase is overwhelming 
in the most of the stellar core. On the whole, in the 
new calculations the ee interactions play an even more 
important role for the physical conditions prevailing in 
the interiors of low-mass RGB stars than had previously 
been suspected (Catelan et al. 1996; Catelan 2005). 

Figure 4 shows the ratio of our present conductive 
opacities to various approximations employed previously. 
One can notice values that are significantly different from 
unity. The difference between our results and Potekhin 
(1999) (solid curves) is due to the neglect of the ee scat- 
tering by the latter. In the case of the comparison with 
Hubbard & Lampe (1969) (dot-dashed curves) the dif- 
ference is mainly due to the deficiency of their treatment 
of strongly coupled and relativistic plasmas, discussed 
above. The comparison with Itoh et al. (1984) (dashed 
curves) shows differences that are partly due to the ee 
collisions, but the most striking difference at the end of 
each curve (that corresponds to the boundary of the He 
core) occurs because T > Tp there, whereas the Itoh et 
al. formulae are valid only for T < Tf. 

3. THE EFFECT ON THE EVOLUTION OF LOW-MASS 
STARS 

In the following we analyze the effect of our new con- 
ductive opacities on the evolution of low-mass stars. For 
the tests involving RGB and HB models we employed the 
same stellar evolution code and input physics discussed 



in Pietrinferni et al. (2004), the only difference being 
the use of our new conductive opacities in place of the 
Potekhin (1999) ones. More in detail, we have computed 
the evolution of the following models: i) M to t = 0.8 M©, 
Z = 10~ 4 , Y = 0.245, ii) Af tot = 0.8 M Q , Z = 10~ 3 , 
Y = 0.246, iii) M tot = 0.85 M Q , Z = 0.004, Y = 0.251, 
and iv) Aftot = 1.0 M , Z = 0.0198, Y = 0.2734. This 
choice of initial stellar masses and chemical composi- 
tion allows us to compare directly these models with the 
results by Pietrinferni et al. (2004), who employed the 
Potekhin opacities. 

3.1. Evolution along the Red Giant Branch 

The evolution of the models has been followed from the 
Zero Age Main Sequence until He ignition at the RGB 
tip. We have included mass loss according to the Reimers 
(1975) formula with n = 0.4. The new conductive opaci- 
ties do not affect the morphology and location of the evo- 
lutionary tracks in the H-R diagram: as expected, only 
the He-core mass at the He-flash and the corresponding 
RGB tip brightness are modified. In Table 1 the main 
properties of the two sets of RGB models computed by 
using alternatively the 'old' k c predictions and the 'new' 
ones are summarized; Fig. 5 displays the differences of 
the model He-core masses and RGB tip brightness as a 
function of Z. 

Our new n c values provide He-core masses at the He- 
flash lower than in the case when the previous k c cal- 
culations are used (model A described in §2.4). The 
difference amounts to 0.006 M@, regardless of the value 
of Z. The reason is that the new k c values are larger 
than the older ones, thus producing a different thermal 
stratification in the He core. The stellar models based on 
the updated k c at the RGB tip are cooler at the center 
of the star but hotter at the point of the maximum off- 
center temperature where the He-ignition occurs. This is 
clearly shown by Fig. 6 (left panel) that displays the ther- 
mal stratification along the He core, at different stages 
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TABLE 1 

Selected properties of the models at the RGB tip. 



version 


^trgb a 


log(L/Lm) b 


logT eff c 


log T r d log p r c 


M cHe f 


M(T max ) s 


loET,„„ h 










0.8 Mq - 


- z = o.oooi - y = 


: 0.245 








new 
old 


12.49 
12.48 


3.2586 
3.2908 


3.6427 
3.6411 


7.863 6.061 
7.870 6.079 


0.4971 
0.5031 


0.178 
0.170 


7.973 
7.968 


-4.085 
-4.164 








0.8 A4© - 


- Z = 0.001 - y = 


0.246 








new 
old 


13.69 
13.68 


3.3384 
3.3696 


3.5958 
3.5946 


7.870 6.020 
7.877 6.039 


0.4852 
0.4913 


0.162 
0.154 


7.969 
7.966 


-4.147 
-4.223 








0.85 Mq 


- Z = 0.004 — Y = 


: 0.251 








new 
old 


13.63 
13.62 


3.3829 
3.4128 


3.5473 
3.5453 


7.873 5.997 
7.880 6.015 


0.4788 
0.4847 


0.167 
0.161 


7.975 
7.970 


-3.933 
-3.989 








1.0 Mq - 


Z = 0.0198 - Y = 


0.2738 








new 
old 


12.52 
12.50 


3.4168 
3.4459 


3.4822 
3.4766 


7.874 5.954 
7.881 5.972 


0.4669 
0.4726 


0.191 
0.128 


7.978 
7.977 


-2.266 
-2.239 





a Age (in Gyr) at the RGB tip. 

b Logarithm of the luminosity (in solar units) at the RGB tip. 
c Logarithm of the effective temperature at the RGB tip. 
d Logarithm of the central temperature (in K) at He ignition. 
"Logarithm of the central density (g cm -3 ) at He ignition. 
f Hc-core mass (in solar units) at He ignition. 

g Location in mass (in solar units) of the off-center temperature maximum at the RGB tip. 
h Logarithm of the maximum off-center temperature (in K) at He ignition. 
'Absolute /-Cousins magnitude of the RGB tip. 



along the RGB, for a selected metallicity. The right panel 
in the same figure shows the correspondence between the 
various internal structures and the location of the models 
along the RGB in the H-R diagram. 

The stellar models based on the new k c achieve the 
thermal conditions required for the He-ignition at a 
slightly higher age (see Table 1) compared to the model 
A mentioned above. This is because the new conductive 
opacities force the stellar structures to have a different 
temperature gradient in the He core. As a consequence, 
the temperature at the base of the H-burning shell is 
slightly lower than for models based on the old k c deter- 
minations (see data shown in Fig. 6). 

The difference in He-core mass at the RGB tip causes a 
change of the RGB tip brightness: new models appear to 
be fainter by A log (I//L Q ) 0.03. This change affects 
distances based on the RGB tip brightness. It is well 
known (see, e.g., Lee et al. 1993; Salaris & Cassisi 1998) 
that the /-band (Cousins) brightness of the RGB tip 

( M TRG B) 

is a powerful distance indicator for old, metal- 
poor and intermediate-mctallicity stellar systems such as 
Local Group Dwarfs, because it is very weakly sensitive 
to the heavy element abundance (and age, for ages above 
a few Gyr). For [M/H] ([M/H] = log(Z/X)-log(Z/X) Q ) 
ranging between —2.0 and —0.6, Mj RGB changes by 
about 0.15 mag. The RGB tip brightness being fixed 
by the He-core mass at the He-flash, any change in M C H e 
translates into a variation in M/ RGB . 

The data in Table 1 show how the new k c decreases 
the /-band brightness of the RGB tip by about 0.07 — 
0.08 mag for metallicities lower than about Z <~ 0.004. 
At solar metallicity, the behavior reverses and the new 
estimate is brighter by <~ 0.03 mag. This reversal of the 
trend is entirely due to the trend of bolometric correc- 
tion BCi with T ff, which overcompensates the smaller 
bolometric luminosity at the RGB tip. In fact, all stellar 
tracks based on the new k c experience the He-ignition at 



fainter bolometric magnitudes and higher effective tem- 
peratures. The higher T e g values are caused by the larger 
(of the order of 0.01 Mq) total mass at the flash, which is 
a consequence of the smaller amount of mass lost along 
the RGB evolution, due to the shorter RGB evolution- 
ary lifetime and lower RGB tip luminosity of the models. 
Because of the trend of BCi with T c g, the final value of 
Mj KGB gets brighter at this high metallicity. 

The M/ RGB values obtained with our new opacities 
turn out to be within about la of the empirical value 
obtained by Bcllazzini, Ferraro, & Pancino (2001), who 
coupled their own detection of the RGB tip in the Galac- 
tic globular cluster ui Centauri (NGC 5139) to the cluster 
distance obtained by Thompson ct al. (2001) from the 
analysis of the cluster eclipsing binary OGLEGC 17. 

Before closing this section, we wish to note that the 
main source of change in the value of M c n e when moving 
from the k c evaluations provided by Potekhin (1999) to 
our new conductive opacities is related to the inclusion 
of the ee interactions in our new treatment of conduc- 
tive transport. More in detail, proper treatment of the 
ee interactions according to the fit (13), but with the 
old treatment of the ee scattering frequency in degener- 
ate matter zA 56 (Urpin & Yakovlev 1980; Potekhin et al. 
1997, 1999), contributes about 70% (i.e., - 0.004 M ) 
of the change in M c h c - The better treatment of i/ e d e c? 
that takes into account the relativistic effect of exchange 
of transverse plasmons (see § 2.2.2) accounts for the re- 
maining one-third of the total variation. 

3.2. Evolutionary Properties during the Core 
He-Burning Phase 

As discussed in the Introduction, any change of the 
He-core mass at He ignition causes a variation in the 
luminosity of the Zero Age Horizontal Branch (ZAHB). 
It is worthwhile to investigate how the change in M c n e 
brought about by the new opacities affects the properties 
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Fig. 6. — Left panel: thermal stratification of the He core at different stages along the RGB evolution of a 0.8 Mq, Z = 0.001, Y = 0.246 
model, computed with the 'old' conductive opacities by Potekhin (1999, dashed lines) and our new computations (solid lines) respectively. 
The hottest sequences correspond to ignition of He at the tip of the RGB. Right panel: location in the H-R diagram of the stellar models 
(labelled with capital letters) whose thermal stratification is shown in the left panel. For each model, the time (in Myr) needed to reach 
the He flash at the RGB tip is also listed. 



of HB stellar models. To this purpose, we computed a 
set of HB models originated from the RGB progenitors 
discussed above. 

We focus our attention on the HB models whose ZAHB 
location is at log(T c ff) = 3.85. We choose this effec- 
tive temperature value° for consistency with our previous 
analysis on this issue. 

The main properties of these models are listed in our 
Table 2. The ZAHB models with the updated 
fainter than those based on the previous estimates, be- 
cause of the lower He-core mass of the progenitor. The 
difference is Alog(L/L Q ) ~ 0.02 within the explored 
mctallicity range. The corresponding V-band ZAHB 
brightness (My (ZAHB)) of the models is reported in Ta- 
ble 2. The new conductive opacities cause an increase in 
M v (ZAHB) ranging from - 0.06 mag at Z = 0.0001 
down to ~ 0.04 mag at solar metallicity. This has the ef- 
fect of slightly decreasing the distance to Galactic globu- 
lar clusters and of systematically increasing their ages by 
^0.7 Gyr. These changes in the globular cluster age and 
distance scales are within current uncertainties on these 
estimates (for more details, see, e.g., Recio-Blanco et al. 
2005; De Angeli et al. 2005, and references therein). 

As to the location in effective temperature of the 
ZAHB models, we have already shown that for a given 
RGB progenitor and for fixed mass-loss efficiency along 
the RGB, the total stellar mass at He ignition increases. 
This occurrence - combined with the decrease of the He 
core mass - will cause the ZAHB model location to be 
slightly cooler with respect to the model based on the 
old opacity predictions. 

A consequence of the lower HB luminosity is that 
the evolutionary lifetime of the core He-burning phase 
is slightly longer, by ~ 5-6%. Together with the 

8 We note that this T e ff value cannot be considered fully rep- 
resentative of the average temperature in the RR Lyrae strip (see 
for instance, Marconi et al. 2003). However, this has no impact on 
the present discussion. 



TABLE 2 
Selected properties of HB models. 



version M 3 . 85 4 log(L/L ) 3 . 85 b M v c 


*HB d 


RGB progenitor: 0.8 Mq - Z = 


= 0.0001 - Y 


= 0.245 


new 0.785 1.734 
old 0.801 1.757 


0.434 
0.377 


83.55 
79.03 


RGB progenitor: 0.8 Mq - Z 


= 0.001 - Y 


= 0.246 


new 0.629 1.645 
old 0.639 1.664 


0.628 
0.582 


94.61 
89.64 


RGB progenitor: 0.85 Mq - Z 


= 0.004 - Y 


= 0.251 


new 0.577 1.575 
old 0.585 1.593 


0.825 
0.750 


102.98 
100.69 


RGB progenitor: 1.0 M - Z = 


0.0198 - Y 


= 0.2738 


new 0.531 1.445 
old 0.538 1.463 


1.066 
1.029 


115.88 
111.53 



a Mass (in solar units) of the ZAHB model located at 
logT eff = 3.85. 

b Logarithm of the luminosity (in solar units) of the 
ZAHB model located at logT eff = 3.85. 

c Absolutc visual magnitude of the ZAHB model lo- 
cated at logT c ff = 3.85. 

d Core He-burning lifetime (in Myr) of the HB structure 
with ZAHB location at logT eff = 3.85. 



lower luminosity of the RGB tip and the fainter ZAHB, 
this might affect the theoretical calibration of the R- 
parameter, i.e., the ratio of the number of HB stars to 
the number of RGB objects brighter than the ZAHB level 
(Iben 1968) that is often used as an indicator of the initial 
He abundance of Galactic globular clusters. Recently, 
Cassisi, Salaris, & Irwin (2003) and Salaris et al. (2004) 
have derived a new calibration of R (using the Potekhin 
1999 conductive opacities) that, applied to a large and 
homogeneous database of empirical values for Galactic 
globulars, provides an initial He-abundance in very good 
agreement with the determination based on recent anal- 
ysis of Cosmic Microwave Background fluctuations and 
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Big Bang nucleosynthesis predictions (see Spergel et al. 
2006 and references therein). Here we have estimated the 
value of the R parameter based on the models presented 
in the previous sections and compared the results with 
ones based on the Potekhin opacities. The use of the new 
k c increases the HB lifetime (<hb) by ~ 5%, but the time 
spent by the RGB progenitor at luminosities higher than 
the ZAHB level also increases by ~ 3.5%. The net effect 
is that, for a fixed initial chemical composition, the value 
of R changes by just +0.03 (~ 2%). The dependence of 
the i?-parameter on the He content is SR/5Y ~ 10 (see 
Cassisi et al. 2003) and therefore the change of the R cal- 
ibration caused by the new k c would lead to a variation 
by only 0.003 of the estimated initial He content of the 
Galactic globular cluster system. This means that the 
reliability of the i?-parameter calibration presented by 
Cassisi et al. and Salaris et al. is not seriously affected 
by the revised conductive opacities. 

3.3. Very Low-Mass Stars and White Dwarfs 

In very low-mass (VLM) stars, i.e. structures with 
mass lower than <~ 0.3 — 0.4 M©, the temperature is of the 
order of the electron Fermi temperature. More in detail, 
considering T/T F w 3.31 x 10- 6 T(/x e /p) 2 / 3 together with 
values of temperature and density characteristic of the 
interiors of VLM stars, one obtains T/Tp <~ 1 — 2. This 
means that the VLM stellar interiors are under condi- 
tions of partial electron degeneracy, the degeneracy level 
strongly decreasing with the total mass. Taking into ac- 
count such circumstantial evidence, we decided to test 
the impact of the new opacity estimates on the struc- 
ture and evolution of VLM stars. We adopted the same 
evolutionary code and physical inputs as in Cassisi et 
al. (2000) but including our new conductive opacities. 
We computed VLM stellar models in the mass range 
0.09 - 0.15 M© for a metallicity Z = 0.001. A close 
comparison with the corresponding evolutionary models 
provided by Cassisi et al. (2000), based on the Itoh et al. 
(1993) conductive opacities, shows that - at least in the 
considered range of stellar mass - there are negligible 
differences between the two sets of stellar models con- 
cerning both the structural properties and evolutionary 
lifetimes. 

The extremely degenerate cores of white dwarfs are 
properly covered by the Itoh et al. and Yakovlev et al. 
opacities, whereas the Hubbard & Lampe (1969) opaci- 
ties are not adequate in this regime (see, e.g., Prada Mo- 
roni & Straniero 2002). However, it is well known (see, 
e.g., Hansen 1999; Salaris et al. 2000; Prada Moroni & 
Straniero 2002) that neither the Itoh et al. nor the Hub- 
bard & Lampe conductive opacities can completely cover 
the He-envelopes of white dwarfs during their entire cool- 
ing phase. A solution generally used (see, e.g., the mod- 
els by Hansen 1999; Salaris et al. 2000) is to employ the 
Itoh et al. (or Yakovlev et al.) opacities, supplemented 
outside their range of validity (e.g., when the matter is 
only weakly degenerate) by the Hubbard & Lampe data. 
Given that the different opacity sources are not based on 
the same assumptions and input physics (see §§1,2) their 
matching in the regime of weak degeneracy is sometimes 
problematic. 

It should be stressed that the Potekhin et al. (1999) 
opacities, calculated a few years ago and aimed espe- 
cially at white dwarf and neutron star applications, have 



already superseded the older results of Itoh et al., as ex- 
plained in §2. Nevertheless, even the Potekhin et al. 
opacities should be corrected following the improvement 
of the treatment of the ee scattering contribution de- 
scribed in §§2.2.2 and 2.3. 

Here we test the impact of the new opacities described 
in § 2 and available at the URL referenced in footnote 7 
on the cooling timescale of white dwarfs. We empha- 
size that the new calculations adequately cover the whole 
white dwarf structure during its cooling sequence. We 
employed the same code and input physics as in Salaris et 
al. (2000) to compute the cooling evolution of a 0.54 M© 
and a 1.0 M© DA white dwarf with a pure carbon core 
(the thickness of the He and H envelope layers is as in 
Salaris et al. 2000) using both our new opacities and the 
combination of Itoh et al. and Hubbard & Lampe (1969) 
opacities mentioned before. The mass-radius relation- 
ship is unchanged, and the cooling times at a given lumi- 
nosity turn out to be almost exactly the same, irrespec- 
tive of the opacities used. The maximum differences are 
of the order of just 1% (the cooling times obtained with 
the new opacities being longer) when log(L/L©) < —5.0 
and the cooling ages are above 14 Gyr. 

4. SUMMARY 

We have reviewed the basic theory to estimate the 
electron-conduction opacity for the stellar matter, and 
discussed limitations and shortcomings of the existing 
calculations that are employed in stellar evolution mod- 
elling. We present new results that can be applied to 
both partial (e.g., in the cores of RGB stars, envelopes 
of white dwarfs) and high degeneracy (e.g., white dwarf 
cores) regimes, for an arbitrary chemical composition. 
Our results update the previous calculations by Hub- 
bard & Lampe (1969), Itoh and coworkers (e.g., Itoh et 
al. 1983, 1984, 1993), and Yakovlev and coworkers (e.g., 
Potekhin et al. 1999; Potekhin 1999; Shternin & Yakovlev 
2006), extensively used in the literature. They improve 
upon Hubbard & Lampe (1969) by including an updated 
treatment of both the ei scattering and the ee scatter- 
ing for strongly coupled and relativistic plasmas, and are 
not restricted to the specific mixtures published by Hub- 
bard & Lampe. Differences with Itoh and coworkers are 
mainly due to their neglect of the ee scattering, inaccu- 
rate treatment of the ei scattering at T ~ 100 — 200, and 
the fact that their results do not extend to the regimes 
of partial degeneracy (T > Tp) and weak ion coupling 
(F < 1). Differences with Potekhin (1999) are essentially 
due to his neglect of the ee scattering. Compared to the 
latest results of Yakovlev's group (Potekhin et al. 1999; 
Shternin & Yakovlev 2006), our new opacities differ es- 
sentially in the extension of the ee scattering results to 
the regime of partial degeneracy. We compared our new 
computations with the aforementioned sets of opacities 
in the temperature-density regime of RGB star cores, 
and found non-negligible differences. On the whole, the 
Hubbard & Lampe (1969) results are the most similar to 
ours in the RGB cores, the maximum differences being 
equal to ~ 50% in the central regions of the cores. 

The effect of using our new opacities for RGB and HB 
calculations instead of the Potekhin (1999) data is to de- 
crease the bolometric luminosity of the RGB tip and the 
ZAHB by Alog(L/L ) - 0.03 and Alog(L/L ) - 0.02, 
respectively. The He-core masses at the RGB tip are 



Updated electron conduction opacities 



11 



decreased by ~ 0.006 M@, the total RGB evolutionary 
timescales are almost unchanged, but the time spent 
along the RGB above the HB increases by ^3 %, and 
the HB evolution is longer by ~ 5%. The calibration 
of the i?-parameter as a function of Y remains almost 
unchanged. 

When our new opacities are applied to white dwarf 
models, the resulting mass-radius relationship and cool- 
ing times are essentially unchanged compared to the re- 
sults obtained employing the combination of Itoh and 
Hubbard & Lampe (1969) opacities that until now was 
needed to cover the whole white dwarf structures. The 



impact on Main Sequence models of very-low-mass stars 
is also negligible. 
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